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We present a finite size spin wave calculation on the Heisenberg antiferromagnet on the triangular 
lattice focusing in particular on the low-energy part of the excitation spectrum. For s — 1/2 the 
good agreement with the exact diagonalization and quantum Monte Carlo results supports the 
Cl^' reliability of the spin wave expansion to describe the low-energy spin excitations of the Heisenberg 

model even in presence of frustration. This indicates that the spin susceptibility of the triangular 
antiferromagnet is very close to the linear spin wave result. 
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I. INTRODUCTION 



The antiferromagnetic Heisenberg model represents the prototype of a quantum many body system where the 
Hamiltonian does not commute with the order parameter -the staggered magnetization-. Therefore, quantum fluctu- 
ations are important and may play a crucial role in the ground state properties of the system. In principle they might 
destroy the long-range order of the classical minimum energy configuration depending, basically, on the dimensionality 
' and the topology of the system. 

In two dimensions and for a triangular lattice, quantum fluctuations are also strongly enhanced by frustration so 
that Fazekas and Andersonl^l argued the possible stabilization of a resonating valence bond (RVB) ground state, i.e., 
. a state which does not break the SU{2^ symmetry of the Hamiltonian. However, after a long period of intensive 
Ph ' theoretical and numerical investigationg3~ll3, no definite conclusion has been settled on the ground state properties 
^ , of this system. In fact -unlike the antiferromagnet on the square lattice where there is a general consensus about 
the ordered nature of the ground state even for s = 1/2- in the frustrated cases the lack of exact analytical results 
is accompanied by difficulties in applying stochastic ruiperical methods, as their reliability is strongly limited by the 
well-known sign problem. A few years ago Bernu et a/J13, with a deep analysis of the exact excitation spectra obtained 
for small clusters with the Lanczos technique, already evidenced an antiferrcpjagnetic 120° ordering in the ground 
state. Only very recently, however, quantum Monte Carlo (QMC) calculations^ have allowed to extrapolate to fairly 
large system sizes giving evidence of a quantum Neel order with the order parameter reduced by about 59% from 
Q\ . its classical value. 

' Assuming as established the long-range Neel order of the ground state, in this paper we address the reliability of 
0^ the spin wave (SW) theory as an analytical tool to treat the Heisenbecg, antiferromagnet on the triangular lattice. To 
this purpose, we generalize to the latter model a previously developedtJ SW theory specialized to finite size systems, 
and we study the low-lying excited states. We will show that the agreement with numerical results obtained by exact 
diagonalization and QMC is very good, thus confirming the effectiveness of SW theory in describing the ground state 
properties of the triangular Heisenberg model. 



Q , II. SPIN WAVE CALCULATION 

o 



Several attempts to generalize SW theosy to finite sizes can be found in the literatureEjTll. Here we will follow the 
method proposed in the mentioned Ref.t3 which allows to deal with finite clusters avoiding the spurious Goldstone 
modes divergences in a straightforward way, and, in particular, without imposing any ad hoc holonomic constraint on 



' the sublattice magnetization. 

Assuming the classical Q — (47r/3,0) magnetic structure lying in the x—y plane and applying the unitary trans- 
formation which defines a spatially varying coordinate system (x'—y'—z') in such a way that the x'-axis coincides on 
each site with the local Neel direction, the transformed Heisenberg Hamiltonian reads: 
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7^ = J ^ [cos (Q • (r, - r,)) {Sf Sf + Sf'sf) 

+ sin (Q . (r, - r,)) (Sf sf - Sf'sf) + Sf Sf] (1) 

where J is the (positive) exchange constant between nearest neighbors, the indices i,j label the points and Vj on 
the A^-site triangular lattice and the quantum spin operators satisfy |Sip = s(s + 1). Then, using Holstein-PrimakofF 
transformation for spin operators to order 1/s, Sf = s — alai, Sf = + Oi): — ~ being a and 

the canonical creation and destruction Bose operators, after some algebra the Fourier transformed Hamiltonian 
results: 
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[^kflkflk + ^-Bk(aJ^al^ + aka_k) (2) 



3Js 

k 

where Ed — — 3Js^iV/2 is the classical ground state energy, Ak = 1 + 7k/2, Bk = —3 7k/2, 
Ik— [cos(fca;) +2 cos (/c2;/2) cos (\/3 ^j,/2)] /3 and k is a vector varying in the first Brillouin zone of the lattice. 
The Hamiltonian Ti-swj can be diagonalized for k 7^ 0,±Q introducing the well-known Bogoliubov transformation, 
flk = UkQ!k + WkCkLki with 

^k + Ek \ / „ ^ / ^k — fk 



Uk = ( 2£k / '^k = -sgn(Bk) ( — ) (3) 



where ek — — is the SW dispersion relation. Such a diagonalization leads to 



T^sw - Ed + —;r- ^ (ek - ^k) 



ek(a|,ak +a^k"-k )• (4) 

k#0,±Q 

The Goldstone modes, instead, cannot be diagonalized with this transformation since they become singular for 
k = and k = ±Q. For infinite systems such modes do not contribute to the integrals in Eq. (^, but in the finite 
case they are important and they must be treated separately. By defining the following Hermitian operators 

Qz = - ao) , (5) 

such that, [Qq,, Q/3] = and [Qa, ^sw] — for a, P — x, y, z, the contribution of the singular modes, TisM, in Eq.(|^) 
can be expressed in the form 

nsM = -SJsAo + 3Js^ [Ql + Ql + Ql] . 



Then, taking into account the fact that to the leading order in 1/s, Qa — y/2/Ns, where 5*" are the components 
of the total spin, TLsm may be also rewritten in the more physical form 

UsM = -iJsAo + 3 [{s-f + {syf + [S'f] , 

which clearly favors a singlet ground state (for an even nuHiber of sites) being ^0 positive definite. This result 
is highly non trivial since we have recovered the Lieb-Mattisll3 property which has not been demonstrated for non 
bipartite lattices. Actually, a similar result is obtained by solving exactly the three Fourier components k = 0,±Q 
of the Heisenberg modelc3; however, our treatment allows to construct a formal expression for the SW ground state 
on finite triangular lattices which keeps the correct singlet behavior. In fact, starting from the usual SW ground 
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N 


36 


48 


108 


oo 


Esw/JN 


-0.5497 


-0.5464 


-0.5410 


-0.5388 


Eqmc/JN 


-0.5581(1) 


-0.5541(1) 


-0.5482(1) 


-0.5458(1) 



TABLE I. Ground state energy for tjh« s = 1/2 Heisenberg antiferromagnet on the triangular lattice, as obtained within the 
finite size SW theory and with a QMCO calculation. For A*' = 36 the exact value of the energy is Eq/ JN — —0.5604. 



state, composed by the 120° classical Neel order plus the zero point quantum fluctuations (i.e, zero Bogoliubov 
quasiparticles) , 

|0)= n "k'exp[i^ataLkll^^> 



k740,±Q 



with |A^) ~ \Sf ~ s), the corresponding singlet wavefunction is obtained by projecting |0) onto the subspace 
5 = 0: 



OO pQO 

\^sw)^ I da dpi e"'^-+"^'^y+'^'^'\0) (6) 

J —oo J — oo J —oo 

and reads \ipsw) e^~°''^°'-'i^^"'°"'°^\0) . In particular the singular modos have no contribution to the ground state 
energy while the computation of the order parameter (defined as in Ref.E3), m^, requires their remotion: 

k7^0,±Q 

Even for s — 1/2, the previous SW calculation predicts a very good quantitativje agreement with exact results on 
small clusters (N < 36) of both ground state energy and sublattice magnptizationE3. For larger lattice sizes (see Table 
H) , a comparison of the SW predictions with the available QMC resultgia also shows a surprisingly good agreement 
for the ground state energy which is conserved also with respect to the extrapolated values in the thermodynamic 
limit. Furthermore the QMC estimate of the order parameter in the thermodynamic limit ^t^qmc ^ 0.41(1) can be 

reasonably compared with the SW prediction — 0.478 . In the next sections we will also show that the low-energy 
spectra on finite sizes, calculated within the SW approximation, compare very accurately with the available numerical 
results. This agreement supports the numerical evidence for an ordered ground state in the present model. 



III. LOW-ENERGY SPIN WAVE SPECTRUM 

In this section, we show how to construct the low-lying energy spectra E(S) for finite systems. Following Ref.0, a 
magnetic field in the z-direction is added to stabilize the desired total spin excitation S, 

T^sw — T^sw — hs E S^. 

i 

Classically, the new solution is the 120° Neel order canted by an angle 9 along the direction of the field h. In order 
to develop a SW calculation, a new rotation around y'-axis is performed on the spin operators and it can be proven 
that takes the same form of Eq.(0) with renormalized coefficients and B]^: 



1 i 2h 2 

2 ^ 2^^' 



Z6J 



being = sm.6. In the present case the only singular mode is k = and its contribution is given by 



SM 



-^Ao^ + 3J§ (5^ 



Ns sin( 



which now favors a value of consistent with the applied field, at the classical level. The Hellmann-Feynman theorem 
relates the latter quantities as it follows: 
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10 20 30 10 20 30 

S(S+1) s(s+i) 

FIG. 1. SW (full dots and continuous line), exact (empty dots and dashed lines) and QMClj (empty squares) low-energy 
spectra as a function of |S^| = S{S + 1) for N = 12, 27, 36, 144 and s = 1/2. 



J_d_ 

Ns dh 



E{h) 



where 



E{h) = 



2h 
'z3J 




..4 



3J£ 

~2~ 



k 



(8) 



and ~ y (^k)^ ~ (^k)^- -'-^ particular, the expansion to order (sh)"^ of the two first terms in Eq.(||) for /i — > 
gives the classical perpendicular susceptibility Xci — 1/9 J, while taking the whole expression the known SW resultEl 
Xsw/Xci = 1 — 0.449/2s is recovered. Finally, in order to evaluate the energy spectrum E{S) of the s = 1/2 case, 
a Legendre transformation E{S) = E(h) + hsS has to be performed. The advantage of this expression is that E{S) 
can be computed within our framework for any size of the lattice and, in particular, in the thermodynamic limit. Of 
course our SW spectra are biased by the underlying classical Neel order which is invariant under the point group C^^ 
and under translations of the magnetic sublattices. 



A. Comparison with numerical results 



The occurrence of symmetry breaking in the ground state for N ^ 00 can be evidenced from the structure of the 
finite size energy spectra. In particular, it is well known that when long-range order is present in the thermodynamic 
limit, the low-lying excited states of energy £'(5') and spin S are predicted to behave as the spectrum of a free quantum 
rotator (or quantum top) 



E{S) - -Bo oc 



5(5*+ 1) 
N 



(9) 



as long as S" <C VAL-jFig. ^ shows £^(5*) vs S{S + 1) calculated within the SW theory compared with the exactt3 (for 
N < 36) and QMCEfl (for 5 = 3, up to = 144) results, being the two last also consistent with the above mentioned 
spatial symmetries of the classical Neel stateE2l. Remarkably SW theory turns out to be efficient to reproduce the 
low-energy spectrum in the whole range of sizes. Furthermore, we can extend our calculation to the thermodynamic 
limit and observe easily the collapse of a macroscopic number of states with different S to the ground state as TV — s- cx). 
This clearly gives rise to a broken SU{2) symmetry ground state, as expected within the SW framework. 
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FIG. 2. Size dependence of 1/2xl (a) and of 1/2xl/2 — l/2xi (b) obtained according to Eq. ( [ic| ) using the (s = 1/2) SW 
excitation spectra. The continuous hne is a quadratic fit for L < 18 in (a) and a guide for the eye in (b). 

B. Spin susceptibilities and anomalous finite size effects 

Whenever the quantum top law is verified, the quantity 

[2xsr^ ^NE{S)[S{S+\)]-' , (10) 

should approach the physical inverse susceptibility l/2xsw for infinite size and for any spin excitation S <^ N . This 
feature is clearly present in the SW theory and it is shown in Fig. |^(a) where the 1/2x5 is plotted ior S = L = \/N 
and approaches the predicted value (l/2xsw = 8.167), even if the correct asymptotic scaling 1/2xl — l/2xsw + 
a/L + b/ L"^ turns out to be satisfied only for very large sizes (L > |3fi). Such feature is also shared by the Heisenberg 
antiferromagnet on the square lattice where a similar SW analysisEj has allowed to account for the anomalous finite 
size spectrum resulting from an accurate QMC calculation. Furthermore, similarly to the latter case, a non-monotonic 
behavior of 1/2xl/2 — 1/2xl (Fig- ||(b)), which should extrapolate to as 1/L according to the quantum top law, 
persists also in presence of the frustration within the SW approximation and is likely to be a genuine feature of the 
Heisenberg model. 

IV. CONCLUDING REMARKS 

In conclusion, we have applied a previously developed finite size spin wave theory to the Heisenberg model on the 
triangular lattice. Comparison for the case s = 1/2 of the low-energy part of the spin wave spectra with the exact 
and the more recent quantum Monte Carlo results reveals a very good agreement. 

The accuracy of our results on finite sizes indicates that spin wave theory is a reliable analytical approximation 
to describe the ground state properties of the present model in the thermodynamic limit iV — > oo. In particular the 
effectiveness of the spin wave theory in reproducing the finite size spectrum strongly suggests that the value of the 
spin susceptibility should be very close to the spin wave prediction. This represents the main finding of this paper. 
Furthermore we have found, to order 1/s, anomalous finite size effects for the spin susceptibilities for s — 1/2, similar 
to the case of the square lattice. 

Finally, this study provides further support to the recent numerical evidences about the existence of long-range 
Neel order in the ground state of the present frustrated model. 
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